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A single axon can generate branches connecting with plenty synaptic targets. Process 
of branching is very important for making connections in central nervous system. The 
interstitial branching along primary axon shaft occurs during nervous system development. 
Growing axon makes pause in its movement and leaves active points behind its terminal. 
A mathematical model is developed to describe and investigate axonal network branching 
process. The model under consideration describes axonal network growth in which the 
concentration of axon guidance molecules manages axon's growth. We model the interstitial 
branching from axon shaft. Numerical simulations show that in the model framework 
axonal networks branch similarly to real neural networks in vitro. 

Keywords: axon; branching model; interstitial; neuron; activity; growth cone 
1. Introduction 

Processes underlying brain self-organization and operation of the cerebral cortex 
is still unclear. Different theoretical approaches describing complex spatio-temporal 
structures and connections were suggested (Andras, 2003; Chauvet, 2002; Chauvet 
and Chauvet, 2002; Hentschel and van Ooyen, 1999; Pelt van et a/., 2005; Segev and 
Ben- Jacob, 2000). These investigations give new understanding of the phenomena 
observed in neural network formation stage. 

Formation of specific nervous system crucially depends on the possibility of cre- 
ation contacts between sets of neurons. For this reason, the main problem of neu- 
robiology becomes revealing mechanisms for recognition of the target neurons by 
axons. 

* Corresponding author. Mail to fgafarov@yandex.ru 
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Evolution of the neuron's connections is important for correct functionality of 
nervous system. Nervous system development is accompanied by axons branching 
out, making contacts with plenty target neurons (Gibson and Ma, 2011). The im- 
portance of understanding mechanisms of temporal and spatial axon's branching 
when investigating neural networks is yet unknown with precision. 

Axonal branches may be characterized by their morphology, complexity and 
branch function which they produce. One shape of branching is arborization. The 
formation of arborization usually takes place in axon's terminal end in the target 
region and leads to formation of arbor branches. Another shape of branching is 
formed by bifurcation of growth cone. As the result two daughter branches appear. 
They move away from each other as was observed in center sensor projections in 
spinal cord. Besides these types of branching there is an interstitial branching which 
usually is formed at some distance from the axon's terminal and grows toward 
the target neuron from main axon. This form of branching is produced by sensory 
collaterals in the spinal cord and by the descending projections from the cortex 
(Gibson and Ma, 2011). 

In the paper we suggest a mathematical model which describes growth of neural 
network taking into account the axons branching process. The model is based on the 
experimental biological data. We use interstitial branching of the axons - the branch 
appears on the axon's shaft behind growth cone (Bastmeyer and O'Leary, 1996; Dent 
et al.^ 1999, 2003; Kalil et al.^ 2000; Soussi-Yanicostas et al.^ 2002; Szebenyi et al.^ 
1998). The interstitial branching happens not by biffurcation but by reorganization 
of the growth cone because after reorganization the growth cone still continues its 
movement. After reorganization of growth cone, on axon's shaft appears a place 
usually called branching point. Here new branch is formed as a response to external 
molecular signals. The main novelty of our approach is consideration of branching 
process and network growth depending on axon guidance molecules (AGM) con- 
centration. The paper is continuation of our previous investigations (Gafarov, 2006; 
Gafarov et a/., 2009). 

2. Axon branching processes 

Axon branching, along with its growth and guidance is tightly controlled and re- 
quires concordant activity of many processes (Gibson and Ma, 2011). 

Evolving growth cones are directed to the corresponding target neurons. Some 
cortex neuron pathways emerge by the evolution of the interstitial branching from 
the axon's body. The interstitial branching appears from the axon shaft sometimes 
very far behind the main growth cone. The branching process is closely connected 
with delay of the growth cone. The duration of delay can take from few hours to few 
days. In this time the expansion of the growth cone takes place. When the cone starts 
moving again, the part of the filopodium and lamellipodium which remains on the 
axon behind its growth cone organizes interstitial branching. In natural conditions 
cerebral cortical axons innervate their target neurons via interstitial branching from 
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axons shaft (Bastmeyer and O'Leary, 1996; Kalil et a/., 2000; Ruthazer et a/., 2010). 

One possible mechanism of interstitial branch formation is that axon shaft it- 
self replies for signals taken from the target which are independent from the main 
growth cone. Another possibility is that the main growth cone recognizes the target 
and demarcates specific domains of axon for future branching. Pauses of the cor- 
tical neurons growth cones movements take place in regions where axon branches 
appear later. In the pausing time the growth cone reorganizes and forms new growth 
cone. The rests of the reorganized cones are left behind on the axon shaft as active 
filopodium and new branches grow from these points (Szebenyi et al.^ 1998). So, the 
question is not only who recognizes "branching signal" but also what type of signal 
is that. Factors released by target neurons can define the place of axon branching. 
Factors from local tissues can provide an instructive cue to promote branch forma- 
tion along the axon. For example, the nerve growth factor (NGF) molecules may 
trigger the growth of the axons interstitial branches growing (Gallo and Letourneau, 
1998). Similar effect was observed with fibroblast growth factor 2 (FGF-2) on pyra- 
midal neuron's axons, but in this case FGF-2 affected region nearby of the growth 
cone(Szebenyi et a/., 2001). In this case the axon branching of separate cortical neu- 
rons takes place without targets and application of FGF-2 just increases branching 
(Kalil et aZ., 2000). The guidance molecule Netrin-1 may stimulate the formation 
of the cortical axon branching (Dent et a/., 2004). The influence of Netrin-1 on the 
branch process formation stimulates space-bounded Ca^^ processes in axon that 
take place in the moment of branching formation (Tang and Kalil, 2005). More- 
over, the local application the Netrin-1 stimulates local Ca^^ transition process 
accompanied by growth of the stimulated branching (Hutchins and Kalil, 2008). 
At last the protein Anosmin commonly known as KALI can assist formation of lo- 
cal interstitial branches from projection neurons in the mammalian olfactory bulb 
(Soussi-Yanicostas et a/., 2002). 

As was shown by fluorescently labeled microtubules and high resolution images, 
the process of branch formation takes place in the growth cone. The growth cone 
distinguishes the future branching points. In the rests left behind the branches can 
appear later (Kalil et a/., 2000). 

Growth of the interstitial branches is accompanied by reorganization of micro- 
tubules (MT) from constrained sets to diffused sets and transformation of long MT 
to short MT. Long-term reorganizations of MT is accompanied by cancellation of 
some axonal processes and increasing and stabilisation of others. The MT move 
inside axonal growth cones developing interstitial branches. During the branch out- 
growth from the axons shaft the MT are reorganized in more plastic form, they may 
go alone out and may be fragmented. Short MT break into developing branches by 
anterograde and retrograde movements (Dent et al.^ 1999). 

Axons branching process also depends on the neuron's activity which may in- 
crease or decrease the branching processes. The investigations made previously with 
cats show that the blocking of the neural activity by tetrodotoxin significantly de- 
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creases the area occupied by axon arbors of thalamocortical fibres (Herrmann and 
Shatz, 1995). The studies of blocking rat's cortical neurons activity with pharma- 
cological inhibitors reveals the decreasing of the branches number only expansion 
of arbors maintained the same level (Uesaka et al, 2005). Therefore, regulation of 
axon activity can enlarge or diminish the size of arbors as well as it can change the 
complexity of structure. 

3. The model of growth and branching 

Direction of growth cones motion is controlled by diffusible chemicals - axon guid- 
ance molecules (AGM) (Goodhill, 1997). Concentration q = ^(rj — r^, t) of AGM at 
the point rj in the moment t depends on the concentration of AGM released by i-th 
neuron at the point r^. Concentration of AGM obeys standard diffusion equation 

- D^AdC, + kc, = J,(r„t). (3.1) 

Here - the diffusion coefficient, k is the degradation coefficient and A^^ is the 
Laplace operator in d-dimensional space. Quantity Ji{rj^t) means some external 
source. Solution of this equation is well-known for a long time and it has the following 
form: 

Ci(r-ri,t) = J Gd{v-v\t)4{v')dv' + j'^dt' j Gd{v - v\t - t')Ji{v\t')dv\ (3.2) 

where c?(r) - initial distribution of concentration and G(i(r, t) - the Green function 
in dimension d\ 

At the initial moment t = there is no AGM, c^(r) = 0, and the source 

Ji{vi,t) = a5^'^\vj-Vi)ji{t), (3.4) 

is localized on z-th neuron. Parameter a describes the amount of AGM released by 
neuron per second. The function ji{t) is perceptivity of z-th neuron at the moment 
t to AGM, that hereafter we will name as activity. We suppose that ji{t) < 1. 
Therefore, we arrive with the following expression for concentration 

c,(r, - r„ t) = a / dt'Gdir^ - r„t - t%{t^). (3.5) 
Jo 

Because of the activity of the neurons has no relaxation form but should depend on 
the surrounded neurons we adopt the following equation for activity, 

= -m + f (jr\t) + J2 ^ikjkit)], (3.6) 

where f[x) — x9(x) and 9(x) is step function. Here jf^^{t) - external source de- 
pending on the time t, ji{t) - activity of i-th neuron, r - relaxation time of activity. 
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The ujik are weights which define influence of k-th neuron on the z-th neuron and 
may take three value, —1,0 and 1. If cjik = 1 then it is excitatory connection and 
vice versa if = — 1 then it is inhibitory connection, = means no influence. 
Because neuron can't make connection with itself then oja — 0. For ji{t) > j^^ the 
(jjik = — 1 and for ji{t) < j^^ the ojik = 1. Here j^^ is a threshold value for activity. 

The z-th axon growth dynamic is described by vector of the axon's terminal 
which obeys to differential equation 

, N 

k=i 

Axon moves if ji > j^^. Parameter A describes axon's sensitivity. At initial moment 
the coordinates of axon's terminal are equal to the coordinates of its neuron and 
activity ji{t) = 0. Also there are no connections between k-th and z-th neurons and 
all (jJik = 0. 

In the framework of model branching process has two stages. At the flrst stage we 
flnd points on the axon's shaft in which the future branch will grow up and we deflne 
direction for growth. To flnd the points we calculate its growth rate, concentration of 
AGM on its terminal and its length. Branching starts if three conditions are fulfllled: 

(i) Velocity of its terminal should be smaller than some threshold value Vg, 

(ii) Concentration of AGM at the terminal should be in range between Cfyiin and 

(iii) Length of axon's shaft should be greater than threshold value Lf^- 




Fig. 1. Axon with branching points. Here and L52 - lengths of the first and the second 
segments of axon between branching points which should exceed threshold value L^^ in order to 
branching process may start. The interval of time t^^ means the time after what the mechanism of 
branch formation starts from branching point. C is concentration axon guidance molecules at the 
moment of creation of branching point. The Cmin and Cmax mean the interval for concentration 
where formation of branch is possible. 

If these conditions are fulfllled the coordinates of this point are stored. New 
branch starts to growing after some time after growth cone reorganization (Dent 
et a/., 2004; Kalil et aZ., 2000) so, we launch timer for interval of time t^^. New 
branch does not grow during this time. After this period of time concentration of 
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AGM is verified at the branching point. If the concentration of AGM belongs interval 
from Cmin up to Cmax (scc Fig. 1), the second stage starts. The axon's branching 
process depends on the activity and the greater it is, the greater is the intensity of 
branching becomes (Uesaka et al.^ 2007, 2005; Ohnami et al, 2008). Plenty branching 
points appear on the axon's shaft during its growth, but only part of them serve as 
the basis for branches. For this reason we consider this process as random process 
(Koene et al.^ 2009). Branching process starts if the following condition is fulfilled: 

> (3-8) 

where P - a random number obtained by generator of the random numbers < P < 
1 and j is the axon's activity. Observations show that the growth cone "touches" 
the space around itself (Aeschlimann, 2000). In analogy of this we calculate AGM 
concentration close to terminal in some points taking into account Eq. (3.5) (see 
Fig. 2(a)). Coordinates of these points are calculated with help of equations 

n sin 0cos(f , Hi = y^^ rt sin 6> sin (/^ , Zi = Za + n cos 6 , (3.9) 

where x^, yi and Zi - coordinates points on the sphere, x^, Ha and Za - coordinates of 
the axon's terminal and the step of angles 9 and is 5°. Next we calculate the AGM 
concentration at the distance from axon's terminal at point and calculate angle 

7 between shaft of axon and line between point and terminal of axon (see Fig. 
2(b)). If this angle belongs interval from 7mm to 7max we add new branch in the set 
of branches. From this set of branches we select single branch with maximum AGM 
concentration. Therefore, we obtain new branch of axon. 




Fig. 2. The procedure of finding point for new branch in 3-dimensional space, (a) Calculation the 
concentration around axon's terminal, (b) Calculation of angle between shaft of axon and the line 
between terminal of axon and the point of observation. Axon is shown with 3 segments, that were 
growing during the time At, the angle 7 and triangle for calculation angle 7. 
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4. Numerical simulations 

In the framework of model we neglect the geometry of neuron's soma. We consider 
it as sphere or circle depending on dimension of problem under consideration. Also 
we suppose that all neurons are identical with respect their form and behavior, the 
positions are fixed. We will call these spheres as neurons. 

Let us consider axonal network which consists of the 12 neurons and then the 
network composed of 18 neurons. The parameters used in the paper are shown in 
the Table 1. 



Table 1. Parameters used for numerical calculations. 



Parameter 


Value 


Amount of AGM per unit second, (Goodhill, 1997; Gundersen and 


a =10-^^ 


Barrett, 1980) 




AGM diffusion coefficient, (Goodhill, 1997) 




Relaxation time of activity, (Vogels et a/., 2005) 


T = Is 


Threshold parameter, (Segev et a/., 2003) 


= 0.51 


Growth's speed of an axon for branching 


^;^ = 5.10-^^ 


Degradation coefficient 


k = 10"^s-^ 


Radius of soma 


r = 5 • 10~^cm 


Minimum concentration of substance for branching 


^mm — J- 


Maximum concentration of substance for branching 


Cmax = 30^ 


Coefficient describing axons sensitivity and motility, (Rosoff et a/.. 


\ — 4 1 


2004) 





In the first experiment we used network from 12 neurons with different distances 
between them. We made experiments 6 times for different values of the external 
activity of neurons jext — 0.05, 0.1, 0.2, 0.3, 0.4, 0.5. At the beginning the activity of 
all neurons was the same and equal jexti which models growth of axonal network 
with blockades such as TTX, APV and DNQX. For each value of the external 
activity we made the numerical experiment 9 times, because the branching process 
has causal origin and connects with neuron activity. We used following parameters: 
initial length of the branch rjj — 1.5-10 cm, angle ^min — 90° and ^max — 

95°, 

the time = 250005 and the threshold for minimal length of axon's segment 
Lfh 5 • 10~^cm. 

In framework of our model neurons may have activity from to 1. Neurons 
with activity around zero are inactive neurons. These neurons both have not been 
stimulated by the external signals and activity of these neurons is depressed by 
chemicals. Therefore, small values of neuron's activity models the axonal network 
in the chemical surround such as TTX, APV and DNQX. 

We calculated amount of branching points and number of formed branches. We 
also obtained concentrations of AGM at branching points, neuron's activity at mo- 
ment of branching, velocity of axon's terminal, length of branch and numbers of 
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neurons between which connections are formed and type of connection. The data 
obtained in 9 numerical experiments are arithmeticahy averaged and represented as 
plots. Except numerical data we show results of numerical analysis axonal network 
with neurons and axons to see the network obtained (see Fig. 6 and 7). 
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Fig. 3. The results of the first experiment. The activity dependence on the number of created 
branches (a), the number of the points of branches (b), and percent of created branches (c). 
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Fig. 4. The results of the second experiment. The activity dependence on the number of created 
branches (a), the number of the points of branches (b), and percent of created branches (c). 
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Fig. 5. The results of the third experiment. The activity dependence on the number of created 
branches (a), the number of the points of branches (b), and percent of created branches (c). 



In the first numerical experiments with activity 0.1 and 0.5 number of branching 
points is 27. This is lowest value except the case with ^ext = 0.05. For activity 
3ext = 0.1 the number of the branches is smaller and for this reason the number of 
the branches is the least. For ^ext — 0.5 activity of neurons is highest. It influences 
velocity of axon's growth and rate of network growth because the concentration of 
AGM depends on the neuron's activity. For great velocity and lower time of growth 
the number of points is smaller. Neurons with activity 2ext = 0.2 have the maximum 
number of branching points, but when activity increases the number of branching 
point decreases (see Fig. 3(b)). Number of generated branches becomes greater with 
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increasing of neuron activity except the case with activity 0.5 (see Fig. 3(a)). It 
is also explained by linear dependence between velocity of growth and activity of 
neurons, the time of growth is smaller than in other cases with smaller activity. 

The ratio of the number generated branches and the number of generated branch 
points is shown in Fig. 3(c). For j^xt = 0.1 only 15% from 27 branching points create 
new branches. The probability of new branch creation from the branching points 
increases up to 67% with increasing of neuron's activity to jext — 0.4. A drastic 
increase of probability to create new branch up to 40% appears for j^xt — 0.2. 

In the second numerical experiment we have 18 irregular located neurons. The 
experiment was made 6 time and results were averaged. Unlike first experiment, we 
used different initial length of branching = 5 • 10~^cm and another disposition of 
neurons. We obtained different results. The number of created branching points and 
branches becomes greater. In Fig. 4(a) we show average number (over 6 experiments) 
of created branches for different neuron's activity. In the first experiment we observed 
drastic increasing of number of created branches at jext = 0.2 with small variations 
for jext > 0.2, but in the second experiment we obtained dependence which looks 
like exponent without jump at point 0.2. 

The number of branching points (see Fig. 4(b)) is distinctive comparing with 
first experiment. In the latter case we observed the maximum at point 0.2 while in 
the second experiment - at points 0.4,0.5. The percent of created branches grows 
exponentially, but for activities 0.3 and 0.4 the results of first and second experiments 
coincides with error in 1% (see Fig. 4(c)). 

Because different results in two experiments were obtained, we repeated first 
experiment but with the length of branching as in the second experiment = 
5 • 10~^cm. We refer this experiment as third experiment. The results are the same 
as in the first experiment (see Fig. 5 and Fig. 3). Therefore, this parameter is of no 
importance for network formation. We found also similar results in all experiments 
namely, (i) the number of created branches for activity jext = 0.1 equals 15% — 16% 
and (ii) for activity jext — 0.3 we obtained 48% — 52% in all experiments. 

5. Discussion 

We developed model to describe the process of axon branching during growth of 
neuron network. The basis of the model is diffusion of AGM, the growth process 
is managed by diffusion equation. We incorporate only one type of branching - 
interstitial branching from the axon's shaft. 

The axon's branching is activity depended process as was shown by (Uesaka 
et a/., 2007, 2005) and (Ohnami et a/., 2008) with help of the pharmaceutical drugs, 
blockades and inhibitors. In the present work we verified that this branching pro- 
cess is close to that obtained in real experiments by (Uesaka et al.^ 2007) who ap- 
plied inhibitors depressing neuron's activity. They observed growing axons and their 
branching process. The inhibitor TTX (tetrodotoxin) depresses the ability of nerve 
fibres to conduct pulses and decreases the neuron's activity. 2.1 zb 0.7 new brunches 
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Fig. 6. Axonal network made in the first numerical experiment with jext = 0.1. (a) The initial state 
of network, (b) The final stage of evolution. There are 11 connections - 8 positive and 3 negative. 
New branches do not form new connections. It is explained in the following way. During the axons 
growth the situation in the axonal environment was not changed: no connections created, the 
gradient of ADM concentration remains the same. For this reason influence on branch terminal and 
axon terminal are the same. On the scale of the initial length of branch the gradient of concentration 
is invariable. 



per one neuron appeared on 14 axons when TTX was applied. In our first numer- 
ical experiment the average number of branches is 2.5 for low activity jext = 0.05 
which models the inhibitor application. The work undertaken by Uesaka et al (2007) 
revealed that 1.5 ±0.4 branches appear when APV/DNQX was applied, that corre- 
sponds to activity lower than jext — 0.05. In case of DNQX inhibitor they obtained 
1.7 lb 0.5 branches which corresponds to 2.25 branches in our numerical experiment 
for activity jext = 0.05. The application of APV in culture gives 5.7 ib 1.2 branches. 
It corresponds activity lying in the range from jext — 0.1 to jext — 

0.2. The results 

of Uesaka et al (2007) without inhibitors correspond the neuron's activity lying in 
the range from 0.3 to 0.4. The second numerical experiment gives the same results 
as the first experiment when TTX, DNQX or APV/DNQX were applied. Without 
inhibitors we obtain neuron's activity lying in the range from 0.1 to 0.2 which differs 
from the first experiment. It may be explained by the difference in the number of 
neurons and their dispositions. Concentration of AGM depends on number of neu- 
rons, distance between them and their activity. The greater the neurons quantity, 
the greater is concentration of AGM. The branching process depends on the AGM 
concentration around the axon's tip and the branch point. 

The expected topology of the axonal network was observed in our model. In 
the framework of the model the AGM represents the single process that participate 
in the axon's branching. In reality the different types of molecules starts different 
branching mechanisms, influence different parts of axon, are observed in different 
parts of nervous system and form specific forms of branching. In our model we do not 
take into account the specific properties of each type of molecules. The branching 
process may be accompanied by contraction of axon branches during growth of 
axonal network and competition of branches at projection to the target axon. In 
these cases there is connection of single branch axon with target neuron, the other 
branches are contracted or change their growth trajectory. 

Realization of these phenomena is not so difficult but the main problem for anal- 
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Fig. 7. Axonal network obtained in the second numerical experiment with jext — 0.5. (a) Initial 
state of the network, (b), (c) Intermediate stages, (d) Final stage. There are 32 connections - 15 
positive and 17 negative. 10 axons made connections without branching, 4 axons made connections 
with help of single branch, 2 axons - with two branches, and 2 axons made connections with help 
of 3 additional branches. 

ysis is time of calculations which becomes very huge even for system made of dozen 
neurons. In future we intend to make calculations for system with hundreds axons 
and take into account different phenomena in the branching process (competitive 
between axonal branches, pruning branches). 
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